Unclassified 


1*.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


24.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No.  0704-0 1 88 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

HDl-TR-2161 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Harry  Diamond  Laboratories 


6c  ADDRESS  (Cfty,  Start,  and  ZIP  Code) 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

U.S.  Army  Laboratory  Command 


Sc  ADORESS  (City,  State,  and  ZIP  Code) 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  Stitt,  and  ZIPCodt) 


8b.  OFFICE  SYMBOL  I  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable!  J 


AMSLC 


to  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NO. 

P6.2.H25  AH25 


1 1 .  TITLE  (Include  Security  Classification) 

The  Potential  and  Electric  Fields  of  a  Conducting  Sphere  in  the  Presence  of  a  Charged  Conducting  Plane 


12.  PERSONAL  AUTHOR(S) 

Clyde  A.  Morrison 


13a.  TYPE  of  REPORT 

Summary 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVEREO 

FROM  Aug  88  TO  Nov  88 


14  DATE  OF  REPORT  (Year,  Month,  Day)  IIS  PAGE  COUNT 

June  1989  I  32 


HDL  project :  1 4291 4,  AMS  code:  61 21 20.H25001 2 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Electrostatics,  spherical  harmonics 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  Identify  by  block  number) 

*  The  solution  is  given  to  Laplace’s  equation  lor  a  conducting  sphere  above  a  ground  plane.  The  solution  includes  the 
presence  of  a  uniform  electric  field  perpendicular  to  the  ground  plane.  The  arbitrary  constants  in  the  solution  are  deter¬ 
mined  by  applying  the  boundary  conditions  on  the  plane  and  on  the  sphere.  Formulas  are  given  for  the  evaluation  of  the 
potential  and  electric  fields  at  an  arbitrary  field  point.  All  the  resulting  formulas  are  expressed  in  a  style  suitable  for  com¬ 
putation.  The  computational  results  will  form  the  basis  of  a  future  comparison  report.  / 

"f  -  U  ■  L-  f 


20.  DISTRIBUTION /AVAILA8ILITY  OF  ABSTRACT 
□  UNCLASSIFIEOrtJNUMITEO  E3  SAME  AS  RPT.  Q  OTIC  USERS 

21  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 

22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Clyde  A.  Morrison 

22b  TELEPHONE  (Include  Area  Code) 
(202)  394-2042 

22c.  OFFICE  SYMBOL 
SLCHD-ST-AP 

OD  Form  1473.  JUN  86 


Previous  editions  are  obsolete. 

1 


Unclassified 


Foreword 


Thirty-some  years  ago  Nick  Karayianis  and  I  tackled  a  problem  in 
electrostatics  which  was  similar  to  the  one  discussed  here.  At  that  time,  I 
was  fresh  from  a  course  in  special  functions  taught  by  E.  D.  Rainville  at 
the  University  of  Michigan  and  I  imparted  my  newly  acquired  wisdom  to 
Nick.  A  few  years  later,  Nick  finished  his  doctorate  at  Indiana  University 
and  returned  the  favor  by  teaching  me  Racah  algebra  and  angular  momen¬ 
tum  theory.  In  that  exchange,  I  received  more  than  I  gave.  During  the 
course  of  “breaking  the  neck”  of  this  problem,  I  frequently  thought  of  our 
earlier  problem  and  missed  Nick’s  able  guidance.  This  report  would  have 
been  much  better  with  his  assistance. 
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1.  Introduction 


In  this  report  we  derive  expressions  for  the  electric  potential  and  electric 
fields  produced  by  a  charged  conducting  sphere  above  a  conducting  plane. 
We  were  unable  to  find  the  problem  done  (potential  and  electric  fields  at 
an  arbitrary  field  point)  in  a  number  of  standard  textbooks  on  electro¬ 
statics  [1-7].*  Further,  textbooks  on  spherical  harmonics  did  not  give  the 
solution  [8-10].  Also,  a  preliminary  search  of  the  literature  revealed  a 
number  of  related  problems,  of  which  two  are  cited  here  [11,  12].  These 
latter  references  approach  the  problem  using  bipolar  coordinates  and  refer 
to  a  number  of  internal  memoranda  which  are  difficult  to  obtain. 

In  most  of  those  references  where  the  geometry  of  the  problem  was  simi¬ 
lar  to  the  one  considered  here,  the  discussion  centered  around  finding  the 
capacitance  of  the  system.  In  a  few  references,  the  problem  of  the  electric 
field  at  a  very  particular  point  (for  example,  the  surface  of  tire  sphere)  was 
considered.  Here  we  are  interested  in  finding  the  potential  and,  conse¬ 
quently,  the  electric  fields  at  an  arbitrary  field  point.  Thus  the  possibility 
of  finding  the  problem  done  in  a  textbook  or  monograph  is  rather  slim,  but 
an  extensive  literature  search  would  quite  possibly  find  the  problem.  Since 
an  extensive  search  would  take  considerable  time,  we  decided  to  do  the 
problem  while  the  seareh  was  underway,  and  as  of  this  date  we  have  not 
located  a  usable  solution.  In  our  analysis,  we  use  spherical  coordinates 
throughout. 


2.  Solution  to  Laplace’s  Equation 

We  consider  a  conducting  sphere  of  radius  a,  at  potential  V ,  located  at  a 
height  h  above  a  conducting  plane.  The  potential  of  the  plane  is  chosen  to 
be  zero.  To  simplify  the  problem,  we  replace  the  conducting  plane  by  an 
identical  sphere  at  potential  -V,  located  at  -h  below  the  plane.  The  coor¬ 
dinate  system  is  shown  in  figure  1. 


* References  are  listed  at  the  end  of  this  report. 
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Figure  1.  Coordinate 
system  of  the  two 
spheres. 


We  consider  the  upper  sphere  first  The  two  linearly  independent  solutions 
to  Laplace’s  equation  (V2<j>  =  0)  in  a  coordinate  system  centered  on  the  up¬ 
per  sphere  are 


r"  P„  (cos  9)  (1) 

and 

Pn  (cos  6)  ,  (2) 


where  n  =  0, 1, 2, ....  and  the  Pn( cos  0)  are  the  Legendre  polynomials  [13]. 
For  convenience,  we  list  the  first  few  Legendre  polynomials: 

P  0(cos  9)  =  1  , 

Px( cos  0)  =  cos  0  , 

P2( cos  0)  =  (3  cos20  -  l)/2  , 

P3( cos  0)  =  [(5  cos2©  -  3)  cos  0]/2  , 
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and  the  general  relations, 


Pn(- cos  6)  =  (-l)"PB(cos  0)  , 

nPn( cos  0)  =  (In  -  1)  cos  0^^  (cos  0)  -  (n  -  l)Pn.2(cos  0)  . 

This  last  relation  is  used  recursively  to  generate  Pn( cos  0)  for  n  >  2  using 
the  n  =  0  and  n  =  1  polynomials  as  initial  values. 

The  solutions  given  in  equation  (1)  diverge  for  large  r  and  are  discarded. 
The  solutions  in  equation  (2)  have  no  singularities  in  the  region  a  <,  r  <  «> 
and  are  retained.  Similarity,  the  retained  solutions  for  the  lower  sphere  are 

^(C°sV)  .  (3) 

The  potential  associated  with  the  uniform  electric  field  is  given  by 


<(>£  =  -Ez 


(4) 


The  general  solution  for  the  potential  is  given  by  multiplying  the  solutions 
given  by  equations  (2)  and  (3)  by  arbitrary  constants  and  summing  over  n 
to  obtain 


4>  =  Z  ^Pn  (cos  0)  +  £  i-P*  (cos  V)  -  Ez  .  (5) 

»=o  r  *=o  A 

The  problem  then  is  to  evaluate  An  and  Bn  by  using  the  boundary 
conditions 


and 


4>*v  , 

for 

r  =  a  , 

4>  =  -v  , 

for 

So 

ii 

& 

o 

II 

for 

o 

II 

N 

(6) 


In  equations  (2)  through  (6)  we  have  avoided  introducing  the  arguments  in 
the  potential  <j>  because  the  particular  variables  will  depend  on  which 
boundary  condition  is  considered. 
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3.  Evaluation  of  the  Constants  An  and  Bn 

The  simplest  boundary  condition  is  the  last  of  equation  (6).  On  the  plane 
z-  0,  we  have  r  =  R  and  0  =  n  -  \jr  so  that 

cos  0  =  -  cos  \j f 


and 

Pn(-  cos  y  =  (-1)"  Pn  (cos  \|/)  • 


Then  from  equation  (5)  we  have 


0  =  Z[(-l)M„+fln] 

H 


Pn (COS  V) 

y/l  +  1  ’ 


(7) 


and  since  each  Pn( cos  y)/r”+1  is  linearly  independent  for  each  n,  we  have 

Bn  =  (-l)"+'An  .  (8) 

To  obtain  explict  values  for  An  we  shall  use  the  first  boundary  condition  in 
equation  (6).  In  order  to  do  this  we  need  the  expansion  of  Pn( cos  \j /)/  Rn+X 
in  terms  of  Pn  (cos  0),  This  expansion  is  given  in  a  number  of  places  [14, 
15]  in  very  general  form,  and  the  result  is  derived  heuristically  in  appen¬ 
dix  A.  Since  we  are  to  use  the  first  boundary  condition  in  equation  (6),  we 
use  the  expansion  given  in  equation  (A- 12)  of  appendix  A  which,  with  ob¬ 
vious  changes,  is 


r" 

(2h )"  +  *  + 1 


P n  (COS  0)  , 


(9) 


where 


is  a  binomial  coefficient. 


Substituting  equation  (9)  into  equation  (5)  gives 

«'■<»  ■  ?  + .?.  *<-»■(’  „+  <c°s  e> 

-  E[h  +  rP\  (cos  0)]  , 
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where  we  have  written 


z-  h  +  r  cos  0 
=  h  +  rP j(cos  0)  . 


From  the  first  of  equation  (6),  the  potential  at  the  surface  of  the  upper 
sphere,  <J)(a,  0)  =  V,  we  obtain 


-.yTjk.  |  V  /  1v.  +  *  +  i  an  + 

»|_a"  +  1  *(  }  (2 hT  +  k  +  \  n  )Ak_ 


Pn (COS  0) 


(ID 


-  E[h  +  aP i  (cos  0)]  , 


where  we  have  used  equation  (8)  to  eliminate  Bk  in  favor  of  Ak.  Since  the 
Pn{ cos  0)  are  orthogonal,  for  each  n  we  can  write 


V  =  U0  -X£(-X)‘{/*  -Eh  , 

0  =  C/,  -  X  L  (-X)*+1  f*  +  -  Ea  , 


and 


(12) 


0  =  Un  -  I(-X)*+"  r  +kk  V  for  n  >  1 


In  equation  (12),  Un  and  X  are  defined  by 


and 


U  = 

Un  a"  +  1 


X  _  A 
2/z 
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The  system  of  equations  given  in  equation  (12)  can  be  written  compactly 
in  matrix  form  as 


F  =  (\-XG)»U  , 


(13) 


where  F  is  the  column  vector  with  components 


Fq  =  V  +  Eh  ,  (14) 

Fx-Ea  , 


(15) 


Fn  =  0 ,  n  >  1  , 

1  is  the  unit  matrix,  with  elements  5^,  G  is  the  symmetric  matrix  with 
components 


and 

The 


Gm  =  (-X)*+f 


n  +m 


(16) 


U  is  the  column  vector  with  components  Un  given  in  equation  (12). 


formal  solution  to  equation  (13)  is 


U  =  B  •  F  , 


(17) 


where  B  is 


B  =  (1  -  \G)-' 

=  1  +  XG+X2G2  +X.3G3  +...  . 


(18) 


The  result  given  in  equation  (18)  can  be  quite  deceptive  in  that  not  all  the 
X  dependence  is  explicit;  the  matrix  G  also  contains  X  as  shown  in  equa¬ 
tion  (16).  All  the  constants  ( V ,  a,  E,  and  h)  are  contained  in  F0  and  Fj,  and 
equation  (17)  determines  each  Un  in  terms  of  these  constants.  Also,  since 
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'  r  T 

B  is  symmetric,  equation  (17)  can  be  written  U  =  F‘  •  B,  if  convenient  (F 
in  the  transpose  vector —  a  row  vector). 

For  a  given  potential,  V,  on  the  sphere,  the  charge  Q  on  the  sphere  can  be 
determined*  from  equation  (12).  By  Gauss'  law,  the  charge  is  determined 
by  A0  or  A0  =  Q  (and  therefore  ( Uo  =  Q/a)).  Hence, 

^  =BqqFo  +BqiF\  , 


and  from  (14), 

Q  =Boo(V +  Eh)+B0,Ea  .  (19) 

Thus,  a  knowledge  of  the  two  components  of  B  determine  the  charge  on 
the  sphere.  We  shall  discuss  equation  (19)  further  when  we  consider  the 
case  when  the  charge  on  the  sphere  is  fixed. 


4.  The  Potential  and  the  Electric  Fields 


Having  determined  the  constants  Un  by  equation  (12),  it  is  necessary  to 
determine  the  potential  <j>  and  the  field  components  Er  (=  -dty/dr),  and 
Fe(=  -3<J)/r38)  using  the  Un.  Substituting  Un  for  An  and  Bn  -  (-1)"  + 
into  equation  (10),  we  have 

r 

u„  - 


^  ^  GnmUn 


Pn-E[h  +  r(P ,)] 


(20) 


£r=I 


+  c/„+„Q  xZc, 


:  U rr 


Pn 


+  EP 1 


(21) 


F9  =X 

n 


un 


’k  S  Gnm  G  r, 


where 


(22) 


*In  MKS  units  replace  Q  by  Qi(4KZ0),  where  4KE0  =  1.112650  X  10  ^  Flm. 
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From  equation  (13)  we  can  write  X  £  G„m  Um  =  Un  -  Fn,  and  using  this  re¬ 
sult  in  equations  (20),  (21),  and  (22)  we  can  express  these  equations  in 
terms  of  Un  and  Fn. 


We  now  let 


(23) 


and 


Yn  = 


n-l-i 


U n  - 
~a  ’">0 


(24) 


Zo  is  set  to  zero  since  To  -0  (see  eq  (27)), 


Zo  =  0  and 


r/.Y 


Z„  = 


«+2 


\n-l~l 
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£/„  n 

-  ,»>  0  . 


(25) 


Using  equations  (23),  (24),  and  (25)  in  (20)  through  (22),  we  have 

<t)(r,  6)  =  X(r)  •  P(0) 

Er(r,Q)=Y(r)>F(Q)  (26) 

Ee(r,d)  =  Z(r)-T(  6) 

for  r  <  2 h. 
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In  equation  (26),  the  components  of  P  are  the  Legendre  polynomials  and 
— ►  — > 

the  components  of  T  are  -dPnfdQ.  The  first  few  values  of  T  are 

To  =0, 

T\  =  sin  0, 

T2  =  3  sin  0  cos  0, 

and  the  recursion  relation  is  ^7) 

T„  -  n(P„-i  -  cos  0  Pn  )/sin  0  . 

The  above  results  are  for  r  <  2 h;  for  r  >  2h  we  have  to  use  the  expansion 
for  Pk( cos  y)/Rk+l  in  terms  of  Pn{ cos  0)  given  in  appendix  A  in  equation 
(A-13).  The  solution  for  Un  given  in  equation  (17)  still  holds,  so  we  can 
write 


<j>  =  £  [Un  -  WH] 

IS 


(28) 


where 


W  -  — - _ 

"  (-X)n  £0 


and  the  field  components  are  given  by 


\n+2 


Er^in  +  DiUn-Wnij 


a 


+  EP 1 


and 


(29) 


£e  =  X  [Un  -  Wn ^  -  ETi  .  (30) 


15 


In  equations  (28),  (29),  and  (30),  we  let 
X6  =  -Eh  , 

X[  =  ^  [(/,  -  ]  -  Er  , 


\«+i 


a 


Xn=\y\  [U„  -W„],n  >  1  , 


(31) 


and 


Y6  =0  , 


* 


Zo=0  , 


-0 


yi+2 

z;=|f)  (c/„-iv„)i,«>i 


Then,  as  in  equation  (26),  we  have 


(32) 


(33) 


<|>  =  *'(r)./»(e)  , 

E t=Y'(r)>P(Q)  , 

E0=Z'(r)*r(0)  , 

forr>2h  . 

— ►  -■> 

The  />(0)  and  T(0)  in  equation  (34)  are  the  same  as  in  equation  (26). 


(34) 
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The  results  given  in  equations  (26)  and  (34)  constitute  a  solution  for  the 
entire  region  a  <  r  <  °°  and  0  <  0  <  7t.  The  procedure  consists  in  finding 
the  coefficients  Un  for  given  V,  a,  h,  and  E  by  using  equation  (17).  After 
this,  we  find  the  components  of  X  (X*),  Y  (F')>  and  Z  (Z')  for  r  <  2h  ( r  > 
2 h)  using  equation  (26)  (eq  (34)). 


As  discussed  following  equation  (19),  when  the  charge  Q  on  the  sphere  is 
fixed,  it  is  necessary  to  determine  V  from  equation  (19)  as 


V- 


Q 

aBoo 


-Eh 


Bo l  F 
BS£a 


(35) 


In  order  to  obtain  V,  we  need  to  evaluate  B ^  and  501,  and  the  result  from 
appendix  B  for  general  B0n  is  given  as 

fion=50„+(-X)"Z^  , 

yAp) 

where 

ap=ap^-x2ap_2 

with 

A0  =  A!  =  1  .  [B-21] 

The  values  of  Bqq  and  B01  can  be  calculated  using  these  results. 

The  solution  to  the  problem  where  the  charge  Q  is  fixed  is  then  obtained 
by  determining  V  by  using  equation  (35).  We  then  insert  this  V  into  equa¬ 
tions  (14)  and  (15)  to  determine  F  and  continue  as  in  the  case  in  which  V 
is  given  (eq  (26)  for  r  <  2 h\  eq  (34)  for  r  >  2 h). 


IB-24] 

[B-20] 


5.  Conclusion 


The  results  derived  here  constitute  a  formal  solution  to  the  potential  and 
the  fields  for  a  conducting  sphere  at  a  given  potential  V  (or  charge  Q)  at 
height  h  above  a  charged  conducting  plane.  Preliminary  calculations  show 
that  the  potential  and  fields  are  sufficiently  represented  by  equations  (26) 
and  (34)  even  for  a  modest  number  of  terms  (—10).  However,  the  region 
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near  r  =  2h  is  not  well  represented  by  these  equations.  The  reason  for  this 
can  be  traced  back  to  equations  (A-12)  and  (A-13).  The  series  represented 
in  these  equations  become  conditionally  convergent  near  r  =  2h  and 
diverge  at  r  =  2 h.  If  it  is  important  to  evaluate  equations  (26)  and  (34)  in 
this  region,  either  analytical  methods  or  computational  methods  will  have 
to  be  devised.  One  way  of  handling  this  problem  is  to  compute,  at  values 
of  8,  a  number  of  values  for  <J>,  Er,  and  Ee  to  be  evaluated  for  r  <2h  and 
r  >  2 h.  An  interpolation  formula  is  used  in  the  immediate  vicinity  of  r  = 
2 h.  The  results  given  in  equations  (26)  and  (34)  will  be  used  in  a  follow¬ 
ing  report  where  we  calculate  4>(r,  0),  Er(r,  6),  and  £e(r,  0)  for  a  number 

of  interesting  cases. 
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Appendix  A.  —  Expansion  of  P.  (cos  \j/)  /  Rk+l  in  terms  of 

P„(cos  0) 

In  the  main  body  of  the  report,  we  used  the  formal  expressions  from  Judd1 
and  Morrison2  for  the  expansion  of  Pk(cos  \\r)/RM  in  terms  of  Pn(c os  0). 
Because  of  the  very  general  nature  of  the  expansion  it  is  useful  to  have  an 
independent  check.  To  do  this,  we  consider  the  expansion  (from  eq  (13)) 

oo  I 

£  =  S  (-1  y  Pe  (COS  0)  r  <x  (A-l) 

'  =  Z  (-1)1  -px  Pe  (cos  0)  r>  x  ,  (A-2) 


where 


R  =  Vjc2  +  r2  +  2xr  cos  0 


(A-3) 


and 


x  +  r  cos  0  =  R  cos  y  . 


(A-4) 


(In  the  final  result  we  let  jc  =  2 h.) 

Our  method  consists  of  taking  successive  derivatives  of  MR  with  respect 
to  x  and  using  (A-4)  to  cast  the  result  into  the  desired  form.  First  consider 


where 


R* 

R 2  ’ 


(A-5) 


'B.  R.  Judd,  Angular  Momentum  Theory  for  Diatomic  Molecules,  Academic  Press,  New  York  (1975). 
2C.  A.  Morrison,  Angular  Momentum  Theory  Applied  to  Interactions  in  Solids,  Springer-Verlag  (1988). 
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From  (A-3)  we  have 


= 


x  +  r  cos  0 


and  using  (A-4)  we  obtain  Rx  =  cos  \j/.  Since P{( cos  y)  =  cos  \|/,  we  have 


d  f 1  (cos \j/) 
dx \R  J-  / ? 


(A-6) 


Repeating  this  process  we  obtain 


2  Pi  (cos  \j/) 
R5 


(A-7) 


d3  A  ^  -  6  P$  (cos  y) 

dx7^?  J=  ¥ 


(A-8) 


In  obtaining  these  results,  it  is  convenient  to  express  all  higher  derivatives 
in  terms  of  Rx: 


Ra  =(1  -R2x)/R  , 

Rm=OR2x-3Rx)IR2  . 

This  procedure  keeps  the  size  of  the  algebraic  expressions  under  control. 

On  the  right  side  in  expression  (A-l)  we  consider  only  the  x-dependent 
part  as 
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1  1 
r~7™ 


^hy(l+m+2)^ 


(A-9) 


^7(i]~-(/+1)(/+2)(/  +  3)^ 


and  consider  the  binomial  coefficient  ^  ^  J  ,  which  is  /,  /  +  1, 

(/  +  1)(/  +  2)12,  (/  +  1)(/  +  2 )(/  +  3)/6  for  k  =  0,  1,  2,  and  3,  respectively. 
So  we  write 


d*  f  1  nt  (cosy) 

d?UJ~(  J)  *■  K*+1  ’ 


(A- 10) 


and 


V 


l  +  k 
l 


l 

,  l+k+\  » 


(A- 11) 


thus,  with  x  set  to  2 h,  we  obtain 


Pk 


■pr —  =  Jj  (-1)'  ^  *  k  j  J7TS7  Pi  (cos  6)  ,  r  <  x.  (A- 12) 


For  *  <  r  in  (A-5)  the  result  corresponding  to  equation  (A- 1 1 )  becomes 


j?*' =*’(,(>*• 


with  k  <  l.  Then  we  obtain,  after  the  substitution  x  =  2 h, 


-  P?  V)  =  |  y  (-ir  fn  P.  (cos  8)  .  r  >  2h.  (A-13) 
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Appendix  B.  —  Fixed  Charge  on  the  Sphere 

When  the  charge  on  the  sphere  is  specified  rather  than  the  voltage,  it  is 
useful  to  write  the  boundary  conditions  given  in  equations  (14)  and  (15) 
as 

U  =  BqoFq  +  Bo\F\ 

U  =  Boo  (V  -  Eh)  +  Boi  Ea  ,  C’} 

and  by  Gauss'  law  the  charge  on  the  sphere,  Q,  determines  A0\  that  is, 

A0  =  Q  .  (B-2) 

Then  from  equation  (12),  UQ  =  QJa  and  from  (B-l)  we  obtain* 

V  =  [Qla  -  (B^h  +  B01a)  E)/Bqo  >  (B-3) 

and  we  notice  that  when  E  =  0  we  obtain 

0  =  cV  ,  (B-4) 


where  c  =  aB0 q. 

As  X  -»  0  {h  -¥  oo)  we  know  that  Bw  =  1  and  the  capacitance  becomes  c  = 
a,  which  is  the  usual  expression  for  the  capacitance  for  a  sphere  of  radius 
a.  From  equation  (B-3)  we  see  that  it  would  be  convenient  to  evaluate  the 
potential  of  the  sphere  and  use  the  results  in  the  calculation  of  the  poten¬ 
tial  Er  and  EQ,  as  is  done  for  the  case  when  V  is  given. 


In  equation  (18)  we  defined  the  matrix  B  as 

B  =  (1  -  XG)"1  ,  (B-5) 


which  we  can  expand  as 


B  =  1  +  XG  +  X2G2  +  -  XpGp  +  •••  ,  (B-6) 


*Ifwe  are  using  MKS  units  we  replace  Qla  by  Ql(4TtZ^i). 
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(B-ll) 


where  y  <  1.  (This  value  of  y  is  much  larger  than  we  need.  In  our  case,  the 
largest  value  of  y  is  1/4.) 
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Using  (B-ll)  in  equation  (B-10)  we  obtain 


(ris_  _  HO" 

G  °"  (i  -X2)"+1  ‘ 


Extending  this  result  to  the  next  term  gives 


(G3)o*  =S(G2)0*G)b, 


-  W  T 

*  (1-X2)* 


and  from  equation  (B-l  1)  we  get 


■3  _  (-X)”(l-X2) 
°"  (1  -X2  -X2)"+5 


(B-12) 


(B-13) 


(B-14) 


We  have  written  the  denominators  in  (B-14)  in  the  particular  form  so  that 
if  we  let  the  denominator  in  (B-13)  be  A2  and  the  denominator  in  (B-14) 
be  Ay  we  have 


For 


a3=a2-\z 


(B-15) 


(G4)*  =I(G3)0*Gjb,  ,  (B-l 6) 


we  have 


(G4  )o«  = 


and  from  equation  (B-ll), 


(B-17) 


(G4)q,  =  (-X)" 


45 


(43  -X242r+1 


(B-l  8) 
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If  we  now  let 


A4  =  A3-X2A2  , 


(B-19) 


we  get  (G*)0n  =  (~^)rt 


A4 


n+T 


The  result  in  (B-19)  suggests  the  form 

,  (B-20) 

with 


A 


o 


=  1 


and 


(B-21) 


Aj  =  l  , 

an  assumption  consistent  with  A2  and  A3  given  above.  To  check  the  result, 
we  calculate  (Gp)0n  using  (B-20).  That  is, 


(G'Jo,  =S(G^1)o 


(Gp)ofl  =  (-XT 


(Vi)" 

(Vi  -^Ap_2r> 


(B-22) 


and  with  the  result  given  in  (B-20)  gives 


0&,  =  W 


(VO" 

(v^ 


(B-23) 
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Then  with  the  result  of  (B-23)  substituted  into  equation  (B-7)  we  obtain 


Bon  =  bon  +  W  X  Xp  ,  (B-24) 

P‘ i  \ap  ) 

and  this  result  with  the  Ap  given  by  the  recursion  relation  in  (B-20)  gives  a 
convenient  form  for  computation  of  Bw  and  £01  for  use  in  equation 
(B-3)  to  determine  the  potential  of  the  sphere. 

It  would  be  nice,  from  an  analytic  point  of  view,  to  be  able  to  calculate 
Bnm  fr°m  (B-6)  by  calculating  each  power  of  G  directly.  However  for 
(GV  we  have 

tVYl 


(C=U=(-V“£x“^^+tm),  QMS) 

but  we  have  not  been  able  to  do  the  sum  in  closed  form.  Karayianis1  has 
shown  that 


Zx2* 

*■ o 


1  - 

(l-Xi)n+m+1 


which  is  Euler’s  transformation  for  the  hypergeometric  function 
F(a,  b\  c ;  X2)  given  by  Rainville2  and  does  not  appear  to  be  capable  of  fur¬ 
ther  simplification. 


‘Nick  Karayianis,  Certain  Summations  Involving  Binomial  Coefficients  and  Their  Relation  to  Dyson's  Conjecture, 
Harry  Diamond  Laboratories,  HDL-TR-1217  (April  1964),  7. 

2E.  D.  Rainville,  Special  Functions,  MacMillan,  New  York  (1960),  Ch  10. 
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